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We report the generation of high-purity twin photon pairs through cavity-enhanced non-degenerate four- 
wave mixing (FWM) in a high-Q silicon microdisk resonator. Twin photon pairs are created within the 
same cavity mode and are consequently expected to be identical in all degrees of freedom. The device is 
able to produce twin photons at telecommunication wavelengths with a pair generation rate as large as 
(3.96 ± 0.03) X 10® pairs/s, within a narrow bandwidth of 0.72 GHz. A coincidence-to-accidental ratio of 
660 ± 62 was measured, the highest value reported to date for twin photon pairs, at a pair generation rate of 
(2.47 ± 0.04) X 10"^ pairs/s. Through careful engineering of the dispersion matching window, we have reduced 
the ratio of photons resulting from degenerate FWM to non-degenerate FWM to less than 0.15. 


High-purity indistinguishable photons underlie many 
quantum optical phenomena^^® and are essential for lin¬ 
ear optical quantum computing.'^’® Their quantum inter¬ 
ference forms the foundation of a vast variety of quantum 
photonic functionalities, such as quantum simulation and 
quantum optical computing, which have been recently 
demonstrated on chip-scale devices.The importance 
of high-purity single-mode indistinguishable photons to 
integrated quantum photonics has excited signihcant in¬ 
terest in the past few years to develop optical fiber and 
chip-scale sources for their generation. 

Here we report the generation of high-purity twin 
photon pairs on a silicon chip. By taking advantage 
of dramatic cavity-enhanced non-degenerate four-wave 
mixing (FWM) in a high-Q silicon microdisk resonator, 
we are able to produce twin photon pairs in a sin¬ 
gle temporal and spatial mode within a very narrow 
bandwidth of 0.72 GHz, with a pair generation rate of 
(3.96 ± 0.03) X 10® pairs/s and high photon-pair cor¬ 
relation, indicated by a coincidence-to-accidental ratio 
(GAR) as large as 660 ± 62 (at a pair generation rate 
of (2.47 ± 0.04) X lO'^ pairs/s), which is significantly 
higher than other approaches reported to date.®®^^® The 
demonstrated device exhibits great potential for on-chip 
quantum photonic applications on the GMOS compatible 
silicon-on-insulator platform. 

Figure 1(a) shows a scanning electron microscopic 
(SEM) image of the device. The silicon microdisk has 
a diameter of approximately 9 pm and a thickness of 
260 nm. A schematic of the pair generation process can 
be seen in Fig. 1(b). Two pump lasers are evanescently 
coupled into the cavity, where twin photon pairs are cre¬ 
ated and then coupled out of the cavity through the same 
tapered fiber. This ensures that the generated photons 
occupy a single spatial mode, which is integral to high 
visibility quantum interference. 
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Figure 2(a) shows the measured cavity transmission 
spectrum of the microresonator. The cavity modes em¬ 
ployed for twin photon generation are quasi-transverse- 
magnetic (quasi-TM) and exhibit high optical Qs ~ 
2.5 X 10® for the loaded cavity modes (Fig. 2(a), in¬ 
sets). The device dispersion was engineered to achieve 



FIG. 1. (a) A scanning electron microscopic (SEM) im¬ 

age of the silicon microdisk sitting on a silica pedestal, (b) 
A schematic illustrating the device, tapered fiber, and non¬ 
degenerate FWM process resulting in twin photon pairs. 
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FIG. 2. (a) A trace of the background normalized cav¬ 

ity transmission. The employed mode family is indicated 
by dashed boxes, with fits of the loaded optical Q shown in 
the corresponding insets. A = Detuning. The calibrated 
free spectral range measurements are depicted between adja¬ 
cent cavity modes, (b) The photoluminescence spectrum as 
a result of simultaneously pumping at the 1496.8 nm and 
1534.3 nm cavity modes, (c) Energy diagrams depicting 
degenerate/non-degenerate FWM. 


an equal spacing among the three modes at 1496.8 nm, 
1514.7 nm, and 1534.3 nm, as indicated by the calibrated 
free spectral range values in Fig. 2(a). With frequency 
matching among the three modes, non-degenerate FWM 
(Fig. 2(c)) can efficiently scatter pairs of pump photons 
from the 1496.8 nm (pump 1) and 1534.3 nm (pump 2) 
cavity modes to create bi-photon states at the central 
mode at 1514.7 nm (Fig. 2(b)), with an efficiency scaling 
as QjV where V is the effective mode volume (see below). 
A simulation using the finite element method (FEM) in¬ 
dicates that the device exhibits a small effective mode 
volume of ^ 10 ^m^. Consequently, the high optical Q 
and small effective mode volume significantly strengthens 
the nonlinear parametric processes inside the cavity, re¬ 
sulting in efficient generation of twin photon pairs. This 


is evident in Fig. 2(b), where the strong photolumines¬ 
cence spectrum clearly demonstrates that photon pairs 
were produced at a single cavity mode with a very clean 
background. 

On the quantum mechanical level, twin photon pairs 
are produced by annihilating two pump photons of dif¬ 
ferent energy eigenstates, resulting in the creation of two 
new photons belonging to the same energy eigenstate. 
This process happens instantaneously inside the cavity, 
in adherence with energy and momentum conservation.^® 
As a result, the created pair of photons are expected to 
be identical in all degrees of freedom such as temporal 
wavepacket, spatial profile, frequency, and polarization, 
which we term as twin photon pairs. This process can 
be understood as the time reversal of degenerate FWM, 
which results in photon pairs with distinguishable fre¬ 
quencies (Fig. 2(c)), and has been employed in numer¬ 
ous works relating to the generation of correlated photon 
pairs in silicon waveguides and resonators. 

Detailed analysis of non-degenerate FWM in a mi¬ 
croresonator shows that the expected probability density 
to emit twin photons from the microresonator into the 
optical fiber, at times ti and t 2 , is given by 

Pciti,t2) « ( 1 ) 

where Fe and Ft are the external photon coupling 
rate and loaded cavity photon decay rate, respectively, 
of the twin-photon mode at frequency ujq- Npi and 
Np 2 are the intracavity photon numbers of the two 
pumps at frequencies ujpi and ujp 2 , respectively, g = 
[cr]n 2 fii^ 0 yyujpiujp 2 ]/(n^V) is the vacuum coupling rate for 
non-degenerate FWM, where rj is the spatial mode over¬ 
lap fraction between the three modes, c is the speed of 
light in vacuum, n 2 is the Kerr nonlinear coefficient, and 
n is the refractive index of silicon. Equation (1) shows 
that the twin photons are strongly correlated in time, a 
clear consequence of the simultaneity of their creation. 
The linewidth of the twin photons is determined by the 
photon decay rate of the loaded cavity. Therefore, the 
twin photons in our device exhibit a narrow linewidth of 
0.72 GHz (corresponding to an optical Q of 2.75 x 10®. 
Fig. 2(a)). The twin-photon generation rate is given by 
Rc = { 8 g‘^NpiNp 2 ){Tl/Tf). Consequently, a high Q/V 
would result in a large twin-photon generation rate and 
long coherence time. 

To characterize the quality of the twin photon pairs, 
we built an experimental setup which is schematically 
depicted in Fig. 3. Two tunable continuous-wave (CW) 
pump lasers were combined in a coarse wavelength divi¬ 
sion multiplexer (CWDM MUX) and evanescently cou¬ 
pled into the silicon microdisk, through an optical ta¬ 
pered fiber. A microscopic image of the device, with 
the tapered fiber fixed on adjacent nanoforks, can be 
seen in the inset of Fig. 3. The nanoforks stabilized the 
transmission signal and allowed for a constant external 
coupling rate for the duration of the experiment. Twin 
photon pairs were generated inside the device and deliv- 
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FIG. 3. A schematic of the experimental setup, used to char¬ 
acterize the twin photon source. MZI: Mach-Zehnder Interfer¬ 
ometer, VOA: Variable Optical Attenuator, CWDM: Course 
Wavelength Division Multiplexing, MUX: Multiplexer, DE¬ 
MUX: De-multiplexer, RPD: Reference Photodetector, LPD: 
Locking Photodetector, TBF: Tunable bandpass filter, and 
SPAD: Single Photon Avalanche Photodiode. The inset shows 
a microscopic image of the device and tapered fiber. 


photon located at the twin-photon mode and its tempo¬ 
rally correlated partner in a separate cavity mode. These 
photons are thus uncorrelated with the twin photon pairs 
and consequently add to the noise background. This is¬ 
sue has presented a major challenge in other schemes for 
twin photon generation.Our device, however, allows 
for the flexible engineering of dispersion which controls 
the phase matching window^^’^® so that it dominantly en¬ 
compasses the desired cavity modes at 1496.8 nm, 1514.7 
nm, and 1534.3 nm (Fig. 2(a)), but degrades considerably 
for the subsequent outer set of cavity modes, thus signif¬ 
icantly reducing the efficiency of degenerate FWM from 
individual pumps. In addition, the tapered-fiber coupling 
scheme enables flexible control of the external coupling of 
cavity modes, which in turn can be applied to control the 
FWM processes involved. As a result, we can optimize 
the taper-disk spacing to maximize the efficiency of the 
non-degenerate FWM process while suppressing the de¬ 
generate ones. Combining these two advantages, we are 
able to produce a highly purified non-degenerate FWM 
process. 

Figure 4 compares the relative strength of the degener¬ 
ate and non-degenerate FWM processes. The degenerate 
FWM contributions were measured in a single pump con¬ 
figuration. The blue photoluminescence spectrum corre¬ 
sponds to pump 1 alone, while the red spectrum corre¬ 
sponds to pump 2 alone, as shown in Fig. 4. To indicate 
their total contribution, the two spectra are numerically 
summed, as seen in the magenta trace. Figure 4 shows 
clearly, that when both pumps are on, the PL spectrum 
(in green) is significantly larger than those produced by 


ered out through the same optical fiber. The transmitted 
composite beam was then passed through a CWDM de¬ 
multiplexer (DEMUX) to separate the twin photon mode 
from the pumps. The twin photon pairs were passed 
through a narrowband (1.2 nm) tunable bandpass filter 
and on to an optical switch, which directed the pairs to 
either a spectrometer or a time-correlated single photon 
counting (TCSPC) setup. The tunable bandpass filter 
was used to cut down noise photons generated from spon¬ 
taneous Raman scattering in the optical fiber.Note 
that twin photons share the same frequency and prop¬ 
agate collinearly, so in order to perform coincidence de¬ 
tection, we passed them through a 50/50 beam splitter 
prior to the coincidence counter. Due to the inherent 
indistinguishability of the twin photons, coincidence de¬ 
tection results in a measurement of the self correlation of 
the bi-photon states. 

One potential issue of twin photon generation is the 
contamination from degenerate FWM processes initi¬ 
ated by individual pump modes. Due to the broad¬ 
band phase matching in micro/nanophotonic devices, the 
non-degenerate FWM process is generally accompanied 
with degenerate processes produced by individual pump 
modes. Consequently, the pump photons may undergo 
degenerate FWM, which produces photon pairs with one 
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FIG. 4. A plot depicting the relative contributions of degen¬ 
erate and non-degenerate FWM to the photoluminescence, at 
a combined dropped optical power of Pd = 37.6 nW. The red 
and blue spectra were generated using a single pump. The 
magenta spectrum results from numerically summing the red 
and blue spectra. The green spectrum was generated with 
both pumps on. 
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FIG. 5. (a) The coincidence spectrum belonging to the last data point in (b). The coincidence counts are in blue with a 

Gaussian fit in red. The accidental coincidence counts are shown in gray with a magnified version depicted in the inset. Note 
that the inset is also plotted as coincidence counts vs. delay time (in ns), (b) A plot of the CAR vs. dropped optical power is 
shown in blue and corresponds to the left vertical axis. A plot of the pair generation rate vs. dropped optical power is shown 
in red and corresponds to the right vertical axis. 


individual pumps. The contribution of degenerate FWM 

can be quantified as ^ ^ ^ ^ , where 

Elpi+p2-iElpi+Elp2) 

Ij {j = pl,p2,pl+p2) is the PL intensity of each process. 
From Fig. 4, we estimate that the ratio of degenerate to 
non-degenerate FWM is '^0.14, clearly showing the dom¬ 
inance of non-degenerate FWM over degenerate FWM. 

The quantum correlation between twin photon pairs 
was characterized by coincidence counting measure¬ 
ments. The photons were detected by two InGaAs sin¬ 
gle photon avalanche photodiodes (SPAD), operated in 
Geiger mode, with a gating frequency of 2.5 MHz, gate 
width of 20 ns, dead time of 10 p,s, and quantum effi¬ 
ciency of 15%. The outputs of the SPADs were trans¬ 
mitted to the TCSPG module for coincidence counting. 
Figure 5(a) shows an example of the coincidence spec¬ 
trum, with a full width at half maximum (FWHM) of 
TpwHM = 840 ps, which is the combined effect of the 
temporal correlation of the twin photons (see Eq. (1)) 
and the instrument response time. The latter was mea¬ 
sured to be about 617 ps, primarily resulting from the 
timing jitters of the two SPADs. To quantify the acciden¬ 
tal coincidences, we shifted the time delay by an integer 
number of the gating period {Tq = 400 ns) and found the 
coincidence spectrum. Figure 5(a) shows an example, in 
gray, as well as a magnified version in the inset. The 
difference between the two spectra clearly indicates the 
high correlation of the twin photons. 

The total coincidence counts Ct were obtained by in¬ 
tegrating the coincidence spectrum over a certain delay¬ 
time window At, with a statistical error given by Ect = 
\/Ct. Accordingly, we can obtain the total accidental 
counts A-t by integrating the accidental coincidence spec¬ 
trum over the same time window. As At is quite small. 


due to the high purity of the twin photons, we first ob¬ 
tained a series of coincidence spectra with a time delay 
of jTo between the two SPADs (j = 1,2, ...,40) and in¬ 
tegrated each of them within the same delay-time win¬ 
dow At to obtain an ensemble of accidental coincidence 
counts Aj. The average and standard deviation of Aj 
gave the accidental coincidence counts At and its statis¬ 
tical error Ea^- Finally, the true coincidence counts were 
obtained by Nc = Ct — At, with a statistical error given 

by En^ = \J(Ect)"^ + {EAt)^■ One important metric 
commonly used to characterize the quantum correlation 
between photon pairs is the coincidence-to-accidental ra¬ 
tio (CAR), which is given by CAR = {Ct — At)IAt, with 

a statistical error of ~ • The 

CAR values and pair generation rate were obtained with 
an integrating window of At = tfwhm, to include the 
major contribution of the coincidence spectrum (Fig. 5). 

The pair generation rate was obtained by calibrating 
the true coincidence counts per gate with the clock fre¬ 
quency, duty cycle, and quantum efficiency of the detec¬ 
tors, as well as the transmittance from the device to the 
detectors for the photon pairs, given as 


, __Ac_ 

“ {mDti)ivt2Dt2)TS’ 


( 2 ) 


where rjti = 0.24 and r]t 2 = 0.25 are the transmittance 
from the microresonator to the SPADs (including the ta¬ 
pered fiber, de-multiplexer, tunable bandpass filter, op¬ 
tical switch and beam splitter, see Fig. 3), Du = 0.15 
is the quantum efficiency of the SPAD, T is the total 
data acquisition time, 6 is the duty cycle of the SPAD 
gating and i = {1,2} refers to the separate twin photon 

















5 


paths. The statistical error of the pair generation rate is 
obtained accordingly from E^c ■ 

Figure 5(b) plots the pair generation rate and CAR 
as a function of the total pump power (including both 
pumps) dropped into the cavity. We recorded a max¬ 
imum CAR of 660 ± 62 with a pair generation rate of 
(2.47±0.04) X 10"^ pairs/s at a dropped power of Pd = 7.0 
fiW. Even with a small combined dropped power of Pd 
= 37.6 we were able to generate twin photons with 
a large pair generation rate of (3.96 ±0.03) x 10® pairs/s 
and CAR of 95 ± 2. Note that these pair generation 
rates were obtained by integrating over tfwhm, which 
omits a significant number of true coincidences due to 
the full width at half maximum of the temporal correla¬ 
tion function, Tc = (21n(2))/ri « 306 ps (see Eq. (1)), 
being comparable to the instrument response time of the 
measurement system, tirf ~ 617 ps. To obtain the to¬ 
tal true coincidences, we increased the integration win¬ 
dow. By varying the size of the coincidence window, 
we found that i?pair saturated at At = 5 x tfwhm, 
which resulted in a maximum pair generation rate of 
(5.29 ± 0.04) X 10® pairs/s. 

In summary, we have demonstrated the generation of 
twin photon pairs in a high-Q silicon microdisk resonator. 
We measured a record CAR of 660 ± 62 and pair gen¬ 
eration rate as large as (3.96 ± 0.03) x 10® pairs/s. The 
twin photon pairs were produced at telecommunication 
wavelengths using micro-Watt power levels on a CMOS 
compatible platform. As the field of quantum photonics 
makes rapid advances toward practical applications, the 
underpinning requirement of high fidelity quantum inter¬ 
ference will remain. 26 jjigh-purity single-mode twin 
photon pairs, which are identical in all degrees of free¬ 
dom and exist within a narrow bandwidth, distinguish 
themselves as a great candidate for being the informa¬ 
tion carriers in such systems. 
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